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We compare predictions of the Capillary Wave Model beyond its 
Gaussian approximation with Monte Carlo results for the energy gap 
and the surface energy of the 3D Ising model in the scaling region. 
Our study reveals that the finite size effects of these quantities are well 
described by the Capillary Wave Model, expanded to two-loop order 
(one order beyond the Gaussian approximation). 
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1 Introduction 

Soft modes play an essential role in the description of finite size effects 
(FSEs) in the fluid interface's free energy (see for instance [|[| and references 
therein). 3D spin systems offer a simple context where these effects appear 
and can be studied, e.g. by using numerical simulations to check theoretical 
predictions. It is well known that, between the roughening and the critical 
temperature of the 3D Ising model, interfaces are dominated by long wave- 
length fluctuations (i.e. they behave as fluid interfaces). Because of these 
fluctuations, a complete control on the description of interfaces, starting 
from the microscopic Hamiltonian, has not yet been reached in the rough 
phase. The usual approach consists in assuming an effective Hamiltonian 
describing the collective degrees of freedom. 

An effective model widely used to describe a rough interface is the cap- 
illary wave model (CWM) |Q. In its simplest formulation one assumes an 
effective Hamiltonian proportional to the variation of the surface's area with 
respect to the classical solution. Because of its non-polynomial nature, un- 
til recently the CWM has been studied only in its quadratic approximation, 
the Hamiltonian being equivalent to a (massless) 2D Gaussian model. The 
CWM has been often identified with the Gaussian model. This model was 
shown to give a good description of several features of the interfaces of the 
3D Ising model, like the logarithmic growth (as a function of the lattice size) 
of the interfacial width in the whole rough phase Q, and the FSEs (as a 
function of the shape of the lattice) of the free interface energy in the scaling 
region of the model Q . 

The corrections that arise when passing from the classical approximation 
to the Gaussian one |5|, ^ y], (which we shall term from now on one-loop 
contributions), only depend on one adimensional parameter (namely on the 
asymmetry u = L2/L1 of the transverse sizes of the lattice) and on the 
boundary conditions. They reduce on symmetric lattices (n = 1) to the well 
known finite-size behaviour of the free energy 0, |^ in the large L limit 

oc o-L^ , 1) 
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a being the reduced interface tension. 

As mentioned above, the CWM has been often identified with its Gaus- 
sian approximation. However, the full CWM is interesting in itself because 
of its simple geometrical meaning (see for instance Ref. y] and references 
therein). Also, it coincides with the Nambu-Goto string action in a partic- 



ular gauge (see the next section). The improvement of Monte Carlo simu- 
lations reached in the last years allows now to test the CWM beyond the 
Gaussian approximation. This is the aim of this paper, in which we investi- 
gate the presence of corrections to the pure Gaussian description of the free 
energy of a rough interface in the scaling region of the 3D Ising model. 

A similar analysis was recently made in the context of the 3D three-state 
Potts model [^] . It was shown that higher order corrections to the Gaussian 
approximation of the CWM (from now on called two-loop contributions) give 
contributions to the functional form of the interface's free energy which can 
be exactly estimated. 

The contributions to the interface free energy (i.e. vacuum diagrams on 
finite volumes) due to higher order expansion of the CWM Hamiltonian are 
finite and can be evaluated in a simple way [|l^]. In this paper we present 
evidence that the results are independent of the regularization used. 

It turns out that the two-loop contributions do not depend only on the 
asymmetry parameter u but also on the adimensional expansion parameter 
proportional to the minimal area of the surface, namely a A = aLiL2. As 
a consequence and in contrast to what happens for the Gaussian model, 
even for symmetric (Li = L2 = L) lattices the finite-size behaviour of 
the free energy defined in Eq. (||) gets corrections proportional to (o"L^)^^ 
which are important in extracting the interface tension value using a fitting 
procedure. In this paper we discuss this picture for the Ising model, checking 
the theoretical prediction by means of Monte Carlo simulations and using 
very different techniques and algorithms. 

Finally, let us remark that the 3D Ising model is related through duality 
to 3D Z2 gauge theory. In particular the physics of interfaces is directly 
linked to the physics of their dual gauge observables, i.e. the Wilson loops, 
and, as a consequence, to the properties of the chromo-electric flux tube in 



the confining phase [12|. This means that all the results that we describe in 
this paper have a direct counterpart in the context of Lattice Gauge Theo- 
ries, and could help to better understand the possible string-like descriptions 
of their infrared behaviour. 

The paper is organized as follows: in Sec. 2 we review the model for rough 
interfaces near the continuum. In Sec. 3 we present an alternative derivation 



of the two-loop calculation given in Ref. |11|. In Sec. 4 we introduce and 
discuss the different observables and Monte Carlo techniques used to test 
the CWM. In Sec. 5 we analyze the numerical results. 



2 The Models 

2.1 The Ising Model 

We consider the 3D Ising model on a regular cubic lattice of size Li, L2 in 
the X-, y-directions, where L2 > Li, and size t in the 2;-direction. In the x- 
and the y-direction periodic boundary conditions are imposed, while in the 
z-direction either periodic or anti-periodic boundary conditions are used, 
depending on the Monte Carlo method that is applied. The Hamiltonian is 
defined by 

H = -J2 ^i^j , (2) 

where the sum is over all nearest-neighbour pairs < ij >, and Si = ±1. The 
corresponding partition function is 

Zi= T. e-^^ , (3) 

{s,=±l} 

where (5 = 1/ {ksT). The critical temperature of the model is estimated 
to be /3c = 0.221652(3) [|l3|. The most precise estimate for the roughening 
couphng is /3^ = 0.4074(3) 0. 

We study the model between the critical and the roughening tempera- 
ture, namely in the region (3r > (5 > (3c- 

While in infinite volume, for P > f3c, the system shows a spontaneous 
symmetry breaking, in finite volume this cannot occur, and interfaces ap- 
pear, separating extended domains of different magnetization. 

In particular, we will consider interfaces that are parallel to the x — y 
plane and fluctuate freely in the third, orthogonal, z-direction. In section 
4 we shall discuss how such interfaces can be generated in a Monte Carlo 
simulation. 

Above the roughening temperature (/?r > /3 > /?c) the step free energy 
of the interface goes to zero. As a consequence the interface behaves es- 
sentially as a 2D critical system: it may be freely translated through the 
medium and the long- wavelength, transverse fluctuations in the interface 
position, i.e. capillary waves, have a small cost in energy (hence cannot be 
neglected in calculations). They can be viewed as the Goldstone modes 
associated with the spontaneous breaking of the transverse translational in- 
variance |jl^. To describe the interface free energy one is therefore forced 
to assume an effective model. On the other side, one has the advantage of 
choosing a Hamiltonian defined directly on the continuum to make analytical 
computations. 



2.2 The Effective Model 

We assume the effective Hamiitonian of tfie interface to be proportional to 
its area. Denoting by Xj(^i, ^2)) ^ = 1,2, 3, the coordinates of a point of the 
interface as functions of the parameters ^a, a = 1,2, with < ^^ < 1, we 
can write the area in the standard reparametrization invariant form, 

A= f' d^i /' dC2 Va , (4) 

JO Jo 

where g = det{gai3), with 

_ dxi dx" 

A system with a Hamiitonian proportional to A coincides with the Nambu- 
Goto model for the bosonic string. The corresponding quantum theory 
is anomalous: depending on the quantization method one finds either the 
breaking of rotational invariance or the appearance of interacting longitudi- 
nal modes (Liouville field). We are, however, interested in interfaces of very 
large size where these difficulties disappear and the rotational invariance is 



restored [!£]; in the infrared limit the theory flows to the massless Gaussian 
model. In order to study the first perturbative correction to this limit it is 
convenient to assume that the main contribution to the interface free energy 
is given in this region by small and smooth deformations of the minimal sur- 
face (which is a flat torus of area L1L2). More precisely, we assume that 
there are no foldings nor self-intersections nor overhangsg. Under these con- 
ditions we can choose as parameters the two longitudinal coordinates x = xi 
and y = X2, by putting ^1 = xi/Li and ^2 = X2/L2; as a consequence, the 
transverse displacement 3:3 of the surface becomes a single-valued function 
of them: X3 = (f){x,y). In such a frame Eq. (^) can be written as 



(6) 




^ From the microscopic point of view this is not obvious a priori: in fact it has been 
observed that the interface at small scales is much more similar to a sponge than to a 
smooth surface |2Cl, fell] . The strong linear correlation between the area and the genus of 
the surface (number of microscopic handles), together with the evidence that the partition 
function summed over all genera behaves like a smooth surface, has led to conjecture that 
a simple non-perturbative renormalization of the physical quantities associated to the 
surface occurs |El|. 



Using these notations, the free energy of a fluid interface can be described 

by 
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(7) 
(8) 
(9) 



where A is an undetermined constant within this approach. The reduced 
Hamiltonian TC is given by 



n[<P]=a {A [0] - L1L2) 



(10) 



In words, TClcj)] is given by the change produced by the deformation (j) in the 
interface's area, measured in units of the interface tension a. g 

One can then take into account the quantum contributions by expanding 
Eq. ( p!o|) in the natural adimensional parameter {aA)~^, A = L1L2. The 
interface tension is the only dimensional parameter of this theory. 

Coming back to the adimensional parameters S,a = Xa/La, and putting 
(p' = \/^(p, the Hamiltonian can be written as 
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where u = L2/L1, and the primes have been omitted. For (crA) ^ 
expanding up to the second order, one obtains 
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(11) 

(12) 

-^ 0, 

(13) 
(14) 
(15) 



•^ Far from the scaling region, the interface tension a in Eq. (y_0|) should be replaced by 
the stiffness K. However, in the scaling region, as the bulk correlation length increases 
and the rotational invariance is restored, K approaches a and for all the values of (3 that 
we studied the difference between them can be safely neglected |3|. 



Retaining only the quadratic term of Eq. (|T^ in the Haniiltonian ([l^, 
one obtains from Eq. @ the one-loop contribution which constitutes the 
Gaussian approximation 

Z\9){u) = ^^{iu)h(i)'^ , (16) 

where rj is the Dedekind eta function 

oo 
v{T)=q'/''Y[{l-qn , g = exp(2^ir) . (17) 

n=l 

This is a well known result in string theory and conformal field theory, and 
coincides^ with the partition function of a 2D conformal invariant free boson 
on a torus of modular parameter t = iu fS], |6|, Q| . 

Within this approximation, the interface partition function takes the 
form 

Z = \ e'""^^^^ Z[3) {u) . (18) 

For u = 1, i.e. Li = L2 = L, one recovers the well known finite size behaviour 
of the interface partition function given in Eq. (|l|) , 



Z = Xe-"^ . (19) 

It has been already verified Q that the inclusion of the one-loop contri- 
bution JT^ ) allows to describe accurately finite size effects on asymmetric 
lattices {u > 1): these turn out to be strong enough to make the classical ap- 
proximation ( p!9[ ) completely inadequate. However, this result does not give 
a definitive answer about the reliability of the CWM hypothesis expressed 
by Eq. (p!o|) : it is in fact well known that the Gaussian model is the fixed 
point of a wide class of possible effective descriptions [^. To identify the 
particular effective Hamiltonian which describes the free fluid interface it is 
then crucial to test higher-order contributions. 



Let us also stress that, among various possibilities, Eq. (10) is the sim- 
plest and most intuitive from a geometrical point of view. Moreover, it does 
not add any new free parameter and, even if this hypothesis is rather old [^, 
it has never been tested beyond the Gaussian approximation until recently 

Pi. 



^The constant i] (i) has been introduced just to normalize Zq (1) 



3 Two-Loop Calculation 

In this section we calculate the contribution to the partition function Z from 
the term given in the Eq. ( p!5| ) which is the first correction to the Gaussian 
Hamiltonian. We write Z in the form 

Z = Aoe-°^Z(9)(n)Zf)(u,aoA^o) (20) 

8aoA Zq •' 



where the expectation value is taken in the free theory. In this section the 
subscript o is added to the bare quantities in order to distinguish them from 
the corresponding renormalized ones. In the Nambu-Goto model Qo = 1. 

Eq. (|2l|) shows that the two-loop contribution can be expressed in terms 
of products of double derivatives of the free Green function. Wick's theorem 
gives 

{j\^if^di2{{y^?)') = [?,u\dlGf + 2,u-\dlGf + 

2dlGdlG + 4(%%G)2]^^^, (22) 



with 



G{z-z') = mme)) , (23) 



where we have introduced the complex variable 2 = ^i + iuS,2 ■ 

Let us denote by {a;} the period lattice, namely the set of points of the 
complex plane of the form {w = m + inu, m, n G Z} . Then the periodic 
boundary conditions can be written simply as 

G{z + uj) = G{z) . (24) 

The Green function should also satisfy 

-A (0(00(0)) = 5(2) (0-1 (25) 

where 



^dei^ude2 



^^ + -7^ ■ (26) 



The — 1 term on the right-hand side of Eq. (^) represents the subtraction of 
the zero mode due to the translational invariance of the interface in the z- 
direction: this is the standard procedure making A invertible in the subspace 



orthogonal to the zero mode (see for instance Ref. [22|). 

The solution of the above equation can be expressed in terms of the 
Weierstrass a function, defined through the infinite product 

a{z) =zT\(\-^ g./..+i{.M2 _ (27) 

We can write the solution as a sum of three terms 

G(z) = -i^ln|a(z)| + !^M!^3^e(z2) + i-(9„^(^))' , (28) 

where E2{t) is the first Eisenstein series 

oo ri 

E2(r) = l-24^-^ , g = exp(27rir) . (29) 

1 1 — a 



The last term in Eq. (p^) accounts for the zero mode subtraction, while the 
first two may be thought of as the real part of an analytic function f{z), 
hence they satisfy A3f?e/(2;) = outside the singularities at z G {u;}. For z 
near a node of the period lattice the second term is regular while the first 
one behaves like — 2;^;: In jz — u;| as it should in order to yield the correct 
normalization of the delta function. The coefficient of the second term is 
uniquely fixed by imposing the periodic boundary conditions of Eq. ( p^ (for 
a different, but equivalent form of the solution see for instance Ref. |2^ p. 
571). 

The double derivative of the Green function can be easily calculated 
using the the formula 

p{z) = -^\n(a{z)) , (30) 

where p{z) is the Weierstrass p-function. We need only the first few terms 
of its Laurent expansion about the origin 

1 vr^ 
p{z) = ^ + —ESu)z^ + . . . , (31) 

where -E'4(t) is the second Eisenstein series. The double pole at the origin 
implies that the limit for ^ ^ ^' in Eq. (p2|) is singular. In order to regularize 



this theory we keep (, — ^' different from zero by putting z 
where e is used as an ultraviolet spatial cut-off. We get 
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with c{u) = {'KuE2{iu)/2> - 1), d{u) = Tr^u^E4{iu)/60 and 
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(32) 



(33) 



The three terms in the square brackets of Eq. ( |32D are cut-off dependent 
quantities. They can be reabsorbed in the renormalization of the couplings 
a,X and g, respectively. In fact, putting a = ao + Sao, A = Aq + 6Xo and 
Q = Qo + Sqo in Eq. (|2^), comparison with Eq. (32) yields 



6ao 
SXo 
6qo 



Qo 



Stt'^£^(To 

\oQoCos{2a)c{u) 

Ane'^ao 
d{u) 



QoCos{Aa 



f{u) 



(34) 
(35) 
(36) 



Of course the renormalized quantities o", A and q should not depend on the 
regularization scheme and in particular on the choice of the parameter a. 
On the other hand, choosing a = ?, we get 5qo = 0, which shows that q 
is not renormalized, at least at the first perturbative order. Thus we can 
safely put the Nambu-Goto value ^ = 1 in the final formula 



A 



'UL1L2 



7] (iu) /rj (i) 



1 + 



fju) 
(TL1L2 



+ 



1 



(aLiL2)2 



(37) 



where /(n) is defined in Eq. (|33|). Note that Eq. (|37| ) is symmetric under 
the exchange Li ^-> L2 as a consequence of the functional relation 



Eo 



T^E2{t) 



6r 



vr 



(38) 



Similarly, it can be shown that also Eqs. (^^36) are symmetric under such 
exchange. Our two-loop result coincides with the one obtained some years 
ago in the context of string theory ||ll|, using the (^-function regularization. 



The fact that two completely independent regularization schemes tell us 
that the coupling q is not renormalized suggests that this a general property, 
even if we have not yet found a rigorous argument to support it. Note that 
in our derivation the splitting between the cut-off dependent part and the 
remainder in Eq. (^) is a crucial point. Such a separation arises in a natural 
way in our regularization, but it is conceivable that different cut-off schemes 
might generate a different splitting. We notice, however, that the functional 
form of the two- loop contribution is preserved by the modular invariance. 
This is not obvious in our derivation because, for sake of simplicity, we dealt 
with a rectangular torus with a purely imaginary modulus r = iu. It is 
possible to develop the theory on a generic torus associated to an arbitrary 
complex modulus r. Any physical quantity must be invariant under the two 
generators of modular group: 

S : T ^ -1/r (39) 

T : T ^T+l. (40) 

The terms in the square brackets of Eq. (|3^ ) are not modular invariant 
because the angle a has not an intrinsic geometric meaning. On the contrary, 
the function F{iu) = f{u) which yields the functional form of the two-loop 
contribution (^) is modular invariant. Its form in a general frame (i.e. r 
arbitrary complex number) is given by 




^-^m{T)E^ (r) - \ 




(41) 



Notice that Eq. (37) has no longer the functional form of the classical ap- 
proximation ( |l9|) even on symmetric lattices: for u = lii can be easily seen 
that Eq. (|^) gives 

where the identity E2{i) = — has been used, which follows directly from 
Eq. (H). 

Let us stress finally that no new free parameter is introduced within this 
approach. 

4 Observables and Monte Carlo Simulations 

Let us discuss the observables of the Ising model that we can use to test the 
functional form of the interface free energy. 



10 



In the finite geometry the degeneracy of the ground state is removed: 
the energy of the symmetric, Z^ invariant, ground state is separated by 
a small energy gap Ai? (or inverse tunneling correlation length) from the 
antisymmetric ground state energy. 

This energy splitting is due to tunneling between the two vacua and 
is directly linked to the free energy of the interface. In the dilute gas 
approximation, in which multi-interface configurations are summed over, 
but interactions between interfaces are neglected, the energy splitting is di- 
rectly proportional to the interface partition function Z, where Z is given 
by Eqs. (|^Jl|,|^. It is then easy to show (see e.g. ^) that 



AE = Z(a,Li,L2) , (43) 

where the usual factor 2 has been reabsorbed into the parameter A appearing 
in the definitions of Z. Let us notice that A has the physical dimensions of 
an energy. In the scaling region of the Ising model all dimensional quantities 
should depend on /3 according to the scaling law 

m) =^ eoo(l-|) ' , (44) 

where ^oo is the bulk correlation length in the continuum limit. The most 
precise estimates for v are in the range from 0.624 to 0.630 [|l^, |2^. Since the 
interface tension is the only dimensional physical quantity appearing in the 
interface free energy, in the following we express all physical observables in 
units of the square root the interface tension ^Ja\^ according to the scaling 
law 

He ' 



a(/3) = V^(^l-^j . (45) 

The measurements of the energy gap Ai?, for different choices of the lattice 
sizes, provides then a first direct check on the functional form of the interface 
free energy and allows one to estimate the interface tension a. To this 
end, we have used two different methods, as explained in the following two 
sections. 

A second observable we have used to check our theoretical prediction is 
the surface energy Eg defined by 

1 3Z 

Es{a,LuL2) = --— . (46) 



11 



This observable is particularly useful because it enables us to isolate explic- 
itly quantum contributions beyond the Gaussian one, as it will be discussed 
below. 

Before describing the method we used to estimate these quantities, let 
us make some general remarks, independent from the observable and MC 
method used, on the range of values of j3 and on the lattice sizes where our 
formulae can be used. 

Applying Eq. (p7|), one should pay attention to avoid spurious effects like 
the incomplete restoration of rotational invariance and the residual presence 
of lattice artifacts, as already discussed. To this end we made our simula- 
tion in the scaling region of the Ising model: the lowest temperature used 
corresponds to /3 = 0.240. (In Tab. 1 we present besides other quantities to 
be introduced below estimates for the bulk correlation length of the Ising 
model at the /3-values used in the present study. The MC estimate is taken 



from Ref. |15]. The other estimates are based on the low temperature series 



that was extended to 15th order by Arisue |16]. We analysed the series with 
the help of inhomogeneous differential approximants |17|. The numbers in 
the table are the results from the [1; 7, 7] approximation.) 

The expansion parameter a A should be small enough to justify the per- 
turbative calculation. 

Interactions between interfaces should be negligible, which means that 
the dilute gas approximation must be satisfied: the dominant contribution 
in the probability of creating an interface is proportional to e~ . Too small 
lattice size can give rise to a high density of interfaces and to non-negligible 
interactions. 

Finally, the smallest size of the lattice (Li in our conventions) should 
not only be greater than the bulk correlation length but also greater than 
the (inverse) deconfinement temperature of the dual gauge model. In fact. 



by duality (see for instance |24], see also [Q]), the broken phase of the Ising 
model on an asymmetric 3D lattice corresponds to the confined phase of the 
3D Z2 gauge model at a finite temperature T^^' = 1/Li. 

Precise information on the finite temperature deconfinement transition 
can be found, for instance, in |2^. (3 is mapped to the gauge coupling 
constant $ = — 2 lntanh/3. For each /? there exists a value Lc{P) such that 
for Li < Lc the gauge system is in the deconfined phase, and Eq. (|37| ) cannot 
be applied. 

The values of Lc for the /3's used are given in Tab. 1. These have been 
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obtained using data taken from Ref. |2(;] and the scaling law 



-^ = T(^HPc-kLc)r , (47) 



where T^^^^ = 2.3(1), Pc ^ 0.7614 and u ~ 0.630 (see |E 

4.1 Energy Gap: the Time— Slice Correlations Method 

The first method we have used to extract the energy splitting AE follows 



the procedure explained in Ref. |^ (see also Refs. |28, |J]); we refer to it as 
the time-slice correlations method (TSC). 

Consider a cylindrical geometry, with t ^ Li,L2 and periodic bound- 
ary conditions in all directions, and define the time-slice magnetization S^, 
where A; = 0, 1, ... , t/2 , as 

^ ^ ni=l n2=l 

with n = (ni, n2, k). If one computes the two-point correlation function 

G{k) = {S^Sk) , (49) 

the low energy levels of the transfer matrix spectrum can be obtained from 
the asymptotic ^-dependence of G{k) 

G{k)-Zj = c2(e-'=^^ + e-(*-'=)^^) + 

c?(e-'=^^' + e-(*-'^)^^')+... (50) 

Zj = l + e-*^^ + ... (51) 

where Zj = tr e~* is the partition function of the Ising model in the transfer 
matrix formalism. 

AE' is the energy of the first (antisymmetric) excited state and turns 
out to be (at least) one order of magnitude greater than AE in the range of 
parameters we have used. The coefficient cq corresponds to the magnetiza- 
tion expectation value and can be used to check the consistency of the MC 
results. 

To perform our MC simulations we used a Swendsen-Wang cluster al- 
gorithm |1. For each value of /? considered, i.e. /? = 0.2246,0.2258 and 
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0.2275, Li and L2 ranged from 10 to 35. We usually fixed t = 120, using, 
for simulations with particularly large Li and L2, bigger sizes t = 240 — 360. 
The values of the energy gap AE, extracted using Eq. (|50[), are reported 
in Tabs. 2, 3 and 4 (for Tab. 4 see also the table caption) Q the confidence 
levels of these fits are always above 70%. 

Using this approach one should pay attention to correlations in MC time. 
In particular it turns out that the two-point correlation functions G{k) 
are affected by very strong cross-correlations in MC time. To take under 
control this problem, we followed the procedure used in Ref. Q, scattering 
the evaluation of time-slice correlations in Monte Carlo time. This has the 
advantage of reducing the cross-correlation matrix to an almost-diagonal 
form and simplify the non-linear fitting procedure. For each /3 and each 
lattice we made about 0.6 — 1.2 • 10^ sweeps (after thermalization) with 
1 — 2 • 10^ measurements/observable. A standard jacknife procedure was 
used to evaluate errors. 

Fig. 1 shows a histogram of the magnetization for a typical lattice size. 
Almost all configurations contain zero or two interfaces, which indicates 
that the dilute gas approximation is respected. The plateau corresponds 
to configurations of the system with two interfaces (because of the periodic 
boundary conditions in the z-direction) . It is easy to see that the presence 
of more interfaces would have the effect of strongly modify it. The two peaks 
correspond to configurations without interfaces. 

4.2 Energy Gap: the Boundary Flip Method 

Another method (which we refer to as the boundary flip (BF) method) 
to evaluate the energy gap AE was introduced by one of the authors in 
Refs. m], ||. 

We consider a system which allows both periodic (p) and antiperiodic 
(a) boundary conditions (be). The partition function of this system is given 
by 

Z = Za + Z, = Y: E ^-^""^'''^^ (52) 

be {Si=±l} 

and the fraction of configurations with antiperiodic boundary conditions is 
given by 

1 = 1 E e-/^^(-) = lE E 4c,ae-^^(^'^'=)=<5,.,>. (53) 

{si=±l} be {si=±l} 



*Some of these data have appeared preliminarily in Ref. 
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An analogous result can be found for perfodic boundary conditions. 

We can express the ratio Za/ Zp as a ratio of observables in this system, 

Za < Olca > /rA\ 

'7~ = ^~R ^ ^^^> 

which turns out to be directly connected to the energy gap l^E. 

To see this, let us express the partition functions of the periodic and 
antiperiodic Ising system in terms of the transfer matrix T. The antiperiodic 
boundary conditions are represented by a spin-flip operator P, which flips 
the sign of all spins in a given z-slice. 

The partition function of the periodic system is given by 

Zp = TrT* , (55) 

while the partition function of the antiperiodic system is given by 

Z„ = TrT*P . (56) 

Since the operators T and P commute, they have a common set of eigenfunc- 
tions. Say the eigenvalues of T are A^ and those of P are pj. The possible 
values oi "Pi are 1 and —1. States that are symmetric in the magnetization 
have Pi = 1 and those that are antisymmetric have pi = —1. The partition 
functions take the form 

Zp = Y.Xi (57) 

i 
and 

Za = Y.Xlpi . (58) 

i 

Let us consider the ratio of the partition functions in the low temperature 
phase. If we assume that 

Aos, Aoa >> Ais , Ala , ••• (59) 

then, ignoring terms of order O (Ais/Aqs) , it is easy to show that 

'^\ = f^ = l-2<4c,a> . (60) 

The interface free energy (inverse of the tunneling mass) is then given by 
AE = -ln(Aoa/Aos) = - ^ In (1 -2 <4c,a>) • (61) 
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Assuming that the number of interfaces is even for periodic boundary condi- 
tions and odd for anti-periodic boundary conditions, the dilute gas approx- 
imation leads to exactly the same relation between the interface free energy 
and the boundary statistics [ pl|] . 

We have used the BF method at /? = 0.2240 for a large number of 
lattices, as reported in Tab. 5. In a few cases we studied two or three values 
of t, in order to check the stability of the results. In general, however, the t 
value needed is smaller than the one we would have needed when using the 
TSC method. For this reason, the BF method is particularly useful near the 
critical point where large Li and L2 must be used. 

For the determination of the ratio of partition functions Za/Zp we em- 
ployed the boundary cluster algorithm proposed by one of the authors 

MM- 

For each simulation given in Tab. 4 (case c) and in Tab. 5 we have made, 
after thermalization, 0.7 — 1.4 • 10^ sweeps, depending on the lattice size. 

Moreover, the BF method bypasses the fitting procedure of Eq. (|50| ) 
required by the TSC method, reducing from t/2 to one (< 6bc,a >) the 
observables needed to evaluate E. This allows to save MC time as well as 
drastically reduces the problems connected to correlations in MC time. 

4.3 Surface Energy 

As stated above, one can also measure the surface energy and compare the 
MC results with the CWM predictions. 

If one assumes Eq. (p7|), it follows from Eq. ( p6| ) {a' = dpa and A' = dpX) 
that 

4^^^^)(a,Li,L2) = a'L,L, - ^ - 4o^/3^f ^ • (62) 



Z, 



Q 



On the other hand, if one uses Eq. (|T9|) or Eq. (18) one gets 

E^^'\a,L^,L2) = a'LiL2-j , (63) 

because one-loop quantum contributions do not depend on (3. In particular, 
from Eq. (37) one expects corrections proportional to {Area)~^ due to the 



two-loop corrections to the CWM which would be absent for a pure Gaussian 
model. 

To measure the surface energy we follow the MC procedure which was 
already used in Ref. |33]. Consider again the ratio Za/Zp. From Eq. ( pO| ) 
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one can express the ratio Za/Zp in terms of the two largest transfer matrix 
eigenvalues 

T^W-¥)^ (64) 

Zjp Z \ AQs / 

where we again ignore terms of order O (Ai^/Aos) • In addition we assume 
that t << 1/AE. This means that for periodic boundary conditions the con- 
figurations without an interface dominate while for anti-periodic boundary 
conditions there is only one interface present in almost all configurations. 
Then, within this approximation, one obtains from Eq. (p^) a definition of 
surface energy which does not depend on the extension t, and 

Es = Ep-Ea , (65) 

with Ep =< H >p and Ea =< H >a. 

We have applied this method at /? = 0.240 for many lattice sizes with t 
ranging from 10 to 40. The surface energies Es for each lattice are given in 
Tab. 6. 

4.4 Local Demon Algorithm 

The variance of the energy stems to a major part from fluctuations on small 
scales. Hence an optimally implemented local algorithm is superior to a clus- 
ter algorithm in solving this particular problem. We used a micro-canonical 



demon algorithm |35|, 36, 37] in combination with a particularly efficient 
canonical update |^8| of the demons. This type of algorithm circumvents 
the frequent use of random numbers. The algorithm is implemented using 



the multi-spin coding technique |p6l 37|. Every bit of a computer word car- 



ries one Ising spin. In order to avoid restrictions of the geometry we simulate 
32 (the number of bits in a word) independent systems. 

We have chosen demons that carry the energies 4, 8, and 16. We take 
one demon for each lattice site. 

First we used the cycle described in [^] to simulate the systems with 
periodic and anti-periodic boundary conditions independently. However, it 
turned out that in contrast to the situation close to the roughening transi- 
tion, near the bulk critical point also frequent updates of the demons carry- 
ing the energies 8 and 16 are needed to obtain autocorrelation times close 
to that of the Metropolis algorithm. As one can see from the CPU times 
summarized in Tab. 1 of Ref. p^, this would mean a considerable increase 
of the CPU time needed. 
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We could, however, overcome this problem by simulating systems with 
periodic and anti-periodic boundary conditions in a single simulation and 
coupling them to the same demon system. It turned out that for our par- 
ticular choice of the update cycle, the integrated autocorrelation time of 
Ep — Ea was almost one order of magnitude smaller than the integrated 
autocorrelation time of Ep and Ea as separate quantities. 

The cycle that we used mostly can be explained as follows. The simula- 
tion was done by performing a cycle of 6 groups, where each group consisted 
of a micro- canonical update of the 32 periodic systems and the demon, the 
32 anti-periodic systems plus the demon and a translation or a shift with 
respect to the 32 copies of the spin system of the demon layer with energy 
4, 8, or 16 (alternating). Each group was finished by updating the demons 
with energy 4. The whole cycle was completed by updating the demons with 
energy 8. We performed a measurement twice in this cycle. 

For the 24 x 24 x 40 lattice at /3 = 0.24 we observed the integrated 
autocorrelation times Tp = 5.3 and r^ = 5.1 of the energy with periodic and 
anti-periodic boundary conditions, while the integrated autocorrelation time 
of the surface energy r^ = 0.6 was about one order of magnitude smaller. 
As unit of the autocorrelation time we took one cycle. The simulation with 
200000 cycles took about 104 h on a SPARC 10 workstation. 

We used the drand48 random number generator from the C-library and 
the G05CAF from the NAGLIB. For some choices of the parameters we 
simulated with both random number generators. The results are consistent 
within the error bars. 

5 Monte Carlo Data Analysis 

In this section we compare the CWM predictions with data extracted from 
MC simulations for the two observables defined in Eqs. (HslE^). 

5.1 FSEs of the Interface Free Energy 

Let us first discuss the fits of the energy gap formula Eq. (p3| ) to the data 
given in Tabs. 2-5, using different functional forms for the interface partition 
function Z, in analogy with Refs. |1^, P9|| . 

We call classical, one-loop and two-loop fits those made with Z given by 



Eq. (|19D, Eq. (|18D and Eq. (|37D, respectively. 

The values of the energy gaps we used range over nearly four orders 
of magnitude, the biggest value being AE = 0.0444 at /? = 0.2246, for a 
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Li = L2 = 13 lattice size, and the smallest /\E = 0.09 • 10~^ at /3 = 0.2275, 
with Li = L2 = 26 (see Tabs. 2 and 4). 

To obtain estimates of AE for large lattice sizes the BF method is more 
efficient than the TSC method. For this reason we have principally used it 
at /? = 0.2240 (Tab. 5), which (from our set of values) is the closest to to the 
critical point. Rather large lattices are needed here (see also Tab. 1). The BF 
algorithm has been also used to take some measurements at /3 = 0.2275, to 
obtain energy gaps for particularly small value of the two-loop contribution 
parameter 

x^f{u)/{aA) , (66) 

where /(u) is defined by Eq. (|33D, with u = L2/L1 and A = -L1L2. 

In Tab. 7 we report the fits we made for each /3: the reduced x^ of 
the classical fits turns out to be always more than one order of magnitude 
greater than those of one and two-loop fits, in agreement with what was 



already shown in |Q, 10 1. Moreover, two-loop fits show a level of confidence 
which is higher than the one-loop fits for all /3's: this is the first evidence 
that the two-loop corrections discussed above correctly describe the data. 

For /5 = 0.2246, 0.2258 and 0.2275 the MC datag are given in Fig. 2, 
while those referring to /5 = 0.2240, which range over a rather different scale 
both in AE and in A, are plotted in Fig. 3. The lines represent the best 
two-loop fits; their values in correspondence of the MC data are also given 
in Tabs. 2-5 {AE^^'^^'^). 

Both the high x^ values of the fits and the clear u-dependence of the 
MC data in Fig. 2 and Fig. 3 show that the classical contribution, Eq. ([T^, 
being only function of A and not of u is completely ruled out. 

Therefore, from now on, we concentrate on the one and two-loop func- 
tional forms. 

Let us now make some further considerations on the different results one 
obtains with the one-loop and the two-loop approximations to the CWM. 
The parameter x, at fixed /3 and for any given lattice, gives the relative 
weight of the two-loop correction with respect to the classical and one-loop 
contributions (see Eq. (|37|)). One expects that a fit on a sample of data 
for which x ^ 1 should give consistent results both at one and two-loop 
and that, as one increases the value of x, the fits made with the two-loop 
contributions should have better confidence levels than those made with the 



^Some data at /3 = 0.2275 and (3 = 0.2240 do not appear in the figures: their behaviour 
is analogous to that of the other data. They have been omitted only in order to have a 
better resolution in the figures. 
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one-loop only. Finally, if the CWM is the correct picture, an upper value of 
X should exists at which higher order corrections should become important 
and fits with the two-loop contribution should start to show higher x^'s- 

The values of x for all samples of data are given, in Tabs. 2-5 in the 
last columns, while the values of a which have been used to evaluate it (see 
Eq. (^)) are given in Tab. 7 (the differences between the one and two-loop 
values are not appreciable in this case). One can see that x ranges from 0.03 
up to the rather large value of 0.27. 

In order to make more stringent fits, we can use all the data we have by 
expressing a through the scaling law (|45|): the free parameters then become 
A/-y/a, o"oo and v. With this assumption Eqs. (^J37) take the form 



Zu = -^V^il-Pc/Pre-^-^^'-^^/f^^'^Z^^Hu) (67) 



a ■ ■ ■ ^ 



^21 



Zu- 



-2v //I Q /m-4:iy' 



^ ^ fiu){l-Pc/P)-"' ^ ^( {I -13c/ (5)- 



CTooA \ (cToo^)" 



• (68) 



According to the usual attitude, we make the fits also using the other scaling 
law 

/3 



ct(/3) = ^^(^-^-Ij (69) 

with the obvious modifications of Eqs. ( p^j6q ). 

We performed the fits on the samples of energy gaps 

{ /^E{x) I x<xcut } , (70) 

starting from the first six data {xcut = 0.058) and adding two data each time, 
so that the last fits contain all the 56 data. The corresponding x^ /dof ob- 
tained with the one-loop (empty circles) and two-loop (full circles) formulae 
written above, as functions of Xcut^ are reported in Fig. 4. 

Prom the figure one can see that the one-loop formula (B^) has a 50% 
confidence level (or more) for x^ut ~ 0.15, which corresponds to 36 data 
over 56 of our sample, while the two-loop formula (^) reaches, with the 
same confidence level, Xcut ~ 0.20, which corresponds to 47 data. The fact 
that both fits on the first sample have a reduced x^ greater than one is 
instead due to the lack of statistics (only 6 data points). For x^ut > 0.2, 
the reduced x^ of the fits made with the two-loop correction indicates that 
higher order corrections could be important. Notice, however, that the 
confidence levels for the two-loop fits are better than those for the one-loop 
fits. The parameter of the fits turns out be very stable, for both types of fits. 
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until the x^ /dof ^ 1: these are given in Tab. 8 for the two hmiting samples 
discussed above. The results obtained for (Too and u are consistent with 
the values reported in the literature and, almost, with each other, the only 
difference being in the constant X/^/a. Let us point out that, if the two-loop 
contributions obtained from the CWM Hamiltonian ([Tl|) are physical, one 



expects the two formulae (67.68|) to give the same constant when the above 
fitting procedure is applied to a sample with Xcut ~ 0. As we can observe 
looking at our sample of data, it is rather difficult to reduce Xcut below 0.05 
without spoiling the reliability of the statistics of our analysis, but we can 
significatively check the Constance of A/^/a making fits on (independent) 
samples in which the x values included are approximately constant. That 
is, we define the sample 

{ AE{x) I Xi<x<x,+, } (71) 

and we denote with Xav the average value of x within each sample. 

The results are given in Tab. 9 and in Fig. 5. It turns out that both the 
ly and a^o values are very stable on the whole sample, discarding the last one 
where Xav is too big. However, while the one-loop and two-loop constants 
are compatible for Xav = 0.04 — 0.07, moving toward greater Xav values 
A/-y/a increases systematically in the one-loop case while remains stable 
when the two-loop formula is used. This indicates that the absence of the 
l/(aA) corrections in the Gaussian formula is artificially compensated by 
the enhancements of its constant, already for Xav ~ 0.1. As a final comment, 
let us note that a comparison with a semiclassical 0^ approach ||8[ suggests 
to fix X/^/a to [4 r(3/4)/r(l/4)] ~ 1.352, in good agreement with the two- 
loop value ~ 1.33(1) of Tab. 8 and Tab. 9. 

5.2 FSEs of the Surface Energy 

In this section we discuss the finite size behaviour of the interface partition 
function from MC measurements of the surface energy defined by Eq. (p^. 
As already said in Sec. 4.4, this approach has the advantage that the 
two-loop contribution appears as an additive correction to the surface free 
energy calculated at the Gaussian (and classical) level, as given by Eq. ( |6^ ) 
and Eq. {pqj, respectively. We decided to make our simulations at /3 = 0.240 
also because very precise estimates of a and a' (two of the three parameters 
on which Eqs. ( p^j6^ ) depends on) are known, 

a = 0.0590(2) , from [3] (72) 

a' = 3.813(2) , from (3|] . (73) 
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Our preliminary step is to assume as external inputs these values and, assum- 
ing they are not biased, to use them to test the reliability of the functional 
form of the second order quantum contribution to the interface free energy. 
Let us define the following surface energy differences (at fixed (3) 

AEs (L) = Es (2L, L/2) - Es {L, L) (74) 

considering two different lattices, one with u = L2/L1 = 4 and the other 
with n = 1, but with the same area A = L1L2 = L?' ■ Then, from Eqs. (63,62) 
one obtains 

AE^g^^L) = (75) 



(CWM) / ^ s _ djjZq 



(20 



AS^"'"'™^(L) 



7(20 



d.zf^ 



7(20 



(76) 



The theoretical prediction of the Gaussian model requires these differences 
to be zero, while that of the CWM at two-loop does not depend any more 
on the third parameter, A'/A. Eq. ([76|) can be written explicitly as 



^Ef^M) ^ ^ [_^ (4) _ ^ (1)] (77) 



O" 



2^2 



where f{u) is given by Eq. (B^) and in the two points we are considering 
takes the values /(I) = 0.25 and /(4) ~ 1.52. 

Using the MC data given in Tab. 6 (taken at /? = 0.240) one can easily 
construct these differences. They are given in the first column of Tab. 10, 
where we have assumed L = 2Li = L2/2. It is clear that all these data are 
not compatible with zero. Assuming the validity of Eqs. (72,73|) and plugging 



the values into Eq. (^) we obtain the theoretical predictions reported in the 
second column of Tab. 10. Monte Carlo data and theoretical predictions are 
plotted in Fig. 6 

A nice Area~^ behaviour is clearly seen in the MC data for L > 16, 
in good agreement with Eq. ([77|). Comparing the last column of Tab. 6 
with Tab. 10 one sees that this corresponds to 2; ~ 0.1, while for x > 0.13 
higher order corrections would be needed and for x > 0.20 the two-loop 
contributions are definitively too small. 

Then we follow the same approach of the preceding section, that is we 
assume different forms of the free energy and choose among them using the 
reduced x^s of the fits with equal number of parameters. In this case, our 
main parameters are fixed from the fits. 
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We fit Eq. (^), which can be written exphcitly as 

Er""^ {L,M) = o'UL, + -^m - \ (78) 

and Eq. (p^). The results are given in Tab. 11, where the parameters refer 
to the two-loop fits while, for the Gaussian fits, only the reduced x^'s are 
given in the second column. 

Finally, we can make an interesting check comparing the ratio of g' and 
o with the equation one obtains from the asymptotic scaling laws. Taking 
the derivative with respect to /? in Eq. (]69|), one gets 



a' 



2a (/3-/3e) • ^^^^ 

In order to compute this ratio we take the last two results of Tab. 11, the 



values of Eqs. (72,73) and assume v ~ 0.63. The results are given in the 
fourth column of Tab. 12. Notice that the RHS of Eq. ( [79|) should also be 
obtained from A'/A: from Tab. 11 one sees that the two estimates agree 
within errors. 

On the other side, one can also assume this ratio as an input to evaluate 
the critical index v, inverting Eq. ([79|) ; one gets 

,^-'(l^-M . (80) 

2cr 

The corresponding values are given in the last column of Tab. 12. The 

agreement with the expected values is good, though the errors are still large. 

6 Conclusions 

In this paper we have demonstrated that the finite size effects of interface 
properties in the 3D Ising model are rather well described by the two-loop 
expansion of the CWM. Our MC data for the energy gap and the surface 
energy prove that there are corrections to the Gaussian approximation, and 
the various fits indicate that these corrections are indeed given (up to even 
higher corrections) by the two-loop result, with no extra coefficients intro- 
duced into the the game. Since there is (yet) no rigorous proof that the 
two-loop corrections are universal (i.e., independent from the regularization 
scheme) we have to consider the predictive power of the two-loop CWM 
model a little bit as a miracle. We consider the present contribution as a 
first step towards a deeper understanding of the physics of rough interfaces 
in terms of effective models. 



23 



References 



[1] 

[2] 

[3] 
[4] 
[5] 
[6] 

[7 



M. P. Gelfand and M. E. Fisher, Physica A 166 (1990) 1 
V. Privman, Int. J. Mod. Phys. C 3 (1992) 857 

F. P. Buff, R. A. Lovett and F. H. Stillinger Jr., Phy.3. Rev. Lett. 15 
(1965) 621 

M. Hasenbusch and K. Finn, Physica A 192 (1993) 342 

M. Caselle, F. Gliozzi and S. Vinti, Phys. Lett. B 302 (1993) 74 

C. Itzykson and J. B. Zuber, Nucl. Phys. B 275[FS17] (1986) 580 

B. Bunk, Int. J. Mod. Phys. C 3 (1992) 889 

M. E. Fisher, J. Phys. Soc. Japan Suppl. 26 (1969) 87 

M. E. Fisher and V. Privman, J. Stat. Phys. 33 (1983) 385 

E. Brezin and J. Zinn-Justin, Nucl. Phys. B 257[FS14] (1985) 867 

G. Miinster, Nucl. Phys. B 340 (1990) 559 

R. K. P. Zia, Nucl. Phys. B 251 [FS13] (1985) 676 

P. Provero and S. Vinti, preprint DFTT 58/93, TPI-MINN 93/46-T, 
|hep-lat 931002^ 



[9 

[lo; 

[11] K. Dietz and T. Filk, Phys. Rev. D 27 (1983) 2944 

[12] M. Casehe, R. Fiore, F. GUozzi, P. Guaita and S. Vinti, preprint DFTT 
57/93, to appear in Nucl. Phys. B 

M. Caselle, F. Ghozzi and S. Vinti, Nucl. Phys. B (Proc. Suppl.) 34 
(1994) 263 

For a review on the subject see F. Gliozzi, Acta Phys. Pol. 23 B (1992) 
971 

[13] C. F. Bailhe, R. Gupta, K. A. Hawick and G. S. Pawley, Phys. Rev. 
B 45 (1992) 10438 

[14] M. Hasenbusch, PhD-thesis, Universitat Kaiserslautern (1992) 

[15] M. Hasenbusch, K. Rummukainen and K. Pinn, in preparation 



24 



[16] H. Arisue, Phys. Lett. B 322 (1994) 224 

[17] M.E. Fisher and H. Au-Yang, J. Phys. A 12 (1979) 1677 

[18] M. Luscher, Nud. Phys. B 180 (1981) 317 

[19] P. Olesen, Phys. Lett. B 160 (1985) 144 

[20] V. S. Dotsenko, G. Harris, E. Marinari, E. Martinec, M. Picco and P. 
Windey, Phys. Rev. Lett. 71 (1993) 811 

[21] M. Caselle, F. Gliozzi and S. Vinti, Nucl. Phys. B (Proc. Suppl.) 34 
(1994) 726 

[22] C. Itzykson and J-M. Drouffe, Statistical field theory, Cambridge 
(1989) 

[23] S. Coleman, Aspects of symmetry, Cambridge (1985) 

[24] J.M. Drouffe and J.B. Zuber, Phys. Rep. C 102 (1983) 1 

[25] J.-C. Le Guillou and J. Zinn-Justin, Phys. Rev. B 21 (1980) 3976 

[26] S. Wansleben and J. Zittarz, Nucl. Phys. B 280 (1987) 108 

[27] S. Klessinger and G. Miinster, Nucl. Phys. B 386 (1992) 701 

[28] K. Jansen, J. Jersak, I. Montvay, G. Miinster, T. Trappenberg and U. 
Wolff, Phys. Lett. B 213 (1988) 203 

[29] M. Caselle, F. Gliozzi, P. Provero and S. Vinti, Nucl. Phys. B (Proc. 
Suppl.) 34 (1994) 720 

[30] R.H. Swendsen and J.-Sh. Wang, Phys. Rev. Lett. 58 86 (1987) 

[31] M. Hasenbusch, J. Phys. (Paris) I 3 753 (1993) 

[32] M. Hasenbusch, Physica A 197 423 (1993) 

[33] M. Hasenbusch and K. Pinn, Physica A 203 189 (1994) 

[34] U. Wolff, Phys. Rev. Lett. 62 361 (1989); Nucl. Phys. B 322 759 
(1989) 

[35] M. Creutz, Phys. Rev. Lett. 50 (1983) 1411 



25 



[36] M. Creutz, G. Bhanot and H. Neuberger, Nud. Phys. B 235 [FSll] 
(1984) 417 

[37] M. Creutz, K.J.M. Moriarty and M. O'Brien, Comput. Phys. Comm. 
42 (1986) 191 

[38] K. Rummukainen, Nud. Phys. B 390 (1993) 621 

[39] N. Ito, Physica A 196 (1993) 591 



26 



p 


Lc 


,- {min) 

^1 


tMC 
?bulk 


tIDA 
?bulk 


0.2240 


12.1 


18 




4.527 


0.2246 


10.5 


13 




3.926 


0.2258 


8.5 


12 




3.170 


0.2275 


6.8 


10 


2.62(2) 


2.556 


0.2400 


3.4 


8 




1.251 



Tab. 1. The first column shows the f3-values used in the present study. 

The corresponding inverse critical temperature Lc (in units of the lattice 

spacing) is given in second column. In the third column appear the 

corresponding minimal sizes used in MC simulations. In the last two 

columns we present estimates for the hulk correlation length from a Monte 

Carlo simulation and from an analysis of the low temperature series with 



inhomogeneous differential approximants jj?|/. 



Li L2 



AE 



^E^CWM) 



13 


13 


0.04437(28) 


0.04441 


0.226 


13 


26 


0.01657(28) 


0.01682 


0.181 


13 


30 


0.01321(57) 


0.01326 


0.196 


13 


34 


0.01069(40) 


0.01059 


0.217 


14 


14 


0.03617(32) 


0.03627 


0.267 


14 


28 


0.01157(62) 


0.01156 


0.156 


14 


32 


0.00916(51) 


0.00880 


0.168 


14 


34 


0.00723(65) 


0.00771 


0.176 


16 


16 


0.02398(29) 


0.02355 


0.149 


18 


18 


0.01451(46) 


0.01468 


0.118 



Tab. 2. The energy gaps AE obtained at (3 = 0.2246 from MC simulations 
with the TSC method are reported together with the best fit values to the 



CWM at 2-loop approximation given by Eq. (jSTj). In the last column the 
values of the two-loop parameter x defined in Eq. (fg^j are given. 
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Li 


L2 


AE 


^^(CWM) 


X 


12 


12 


0.03750(23) 


0.03765 


0.184 


12 


24 


0.01126(13) 


0.01117 


0.212 


12 


26 


0.00939(13) 


0.00940 


0.154 


12 


28 


0.00769(29) 


0.00793 


0.161 


12 


30 


0.00665(21) 


0.00672 


0.170 


13 


26 


0.00662(23) 


0.00684 


0.126 


13 


28 


0.00554(32) 


0.00562 


0.130 


14 


14 


0.02225(16) 


0.02212 


0.135 


16 


16 


0.01236(45) 


0.01222 


0.104 


18 


18 


0.00617(48) 


0.00631 


0.082 



Tab. 3. The same as Tab. 2, but for P = 0.2258. 



Li 


L2 


AE 


^^(CWM) 


X 


10 


10 


0.04334(8)^^ 


0.04335 


0.265 


10 


18 


0.01439(26)^ 


0.01457 


0.131 


10 


20 


0.01099(23)^ 


0.01148 


0.136 


10 


23 


0.00797(23)'' 


0.00813 


0.147 


10 


26 


0.00584(23)'' 


0.00583 


0.162 


10 


28 


0.00396(44)'' 


0.00469 


0.174 


10 


32 


0.00289(32)'' 


0.00307 


0.200 


11 


30 


0.00169(35) 


0.00213 


0.140 


11 


35 


0.00114(38) 


0.00113 


0.164 


12 


12 


0.0217(1)'^ 


0.0217 


0.118 


12 


15 


0.01296(21)'= 


0.01274 


0.099 


14 


14 


0.00989(7)'^ 


0.00979 


0.087 


14 


18 


0.00431(7)^= 


0.00432 


0.071 


16 


16 


0.00400(7)" 


0.00397 


0.066 


18 


18 


0.00151(8)'^ 


0.00144 


0.052 


20 


20 


0.000454(11)'= 


0.000466 


0.042 


24 


24 


0.000035(2)^ 


0.000034 


0.029 


26 


26 


0.000009(1)'^ 


0.000008 


0.025 



Tab. 4. The same as Tab. 2, but for f3 = 0.2275. The (a) and (b) are data 
taken from Ref. j^^ and j^, respectively; the (c) data have been obtained 

with the BF method. 
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Li Lo 



AE 



~Xe^cWmY 



18 


18 


54 


0.0230(2)" 


0.0229 


0.161 


18 


30 


60 


0.008628(95) 


0.008650 


0.124 


18 


35 


70 


0.005993(74) 


0.006046 


0.128 


18 


40 


120 


0.004242(42) 


0.004289 


0.137 


18 


45 


90 


0.003034(41) 


0.003076 


0.149 


18 


50 


100 


0.002214(35) 


0.002224 


0.164 


20 


30 


80 


0.006156(90) 


0.006149 


0.101 


20 


35 


70 


0.004010(72) 


0.004021 


0.101 


20 


40 


80 


0.002681(53) 


0.002669 


0.105 


20 


45 


90 


0.001810(39) 


0.001790 


0.112 


20 


50 


150 


0.001205(14) 


0.001210 


0.121 


24 


24 


72 


0.00657(8)*^ 


0.00646 


0.091 


24 


35 


70 


0.001925(32) 


0.001885 


0.071 


24 


40 


80 


0.001102(22) 


0.001107 


0.070 


24 


45 


180 


0.000652(11) 


0.000657 


0.071 


24 


50 


100 


0.000407(10) 


0.000393 


0.074 


24 


60 


120 


0.000145(5) 


0.000143 


0.084 


30 


30 


90 


0.00136(3)^* 


0.00133 


0.058 



Tab. 5. Data at [3 = 0.2240 obtained with the BF method, where the data 
with index (a) have been taken from Ref. jgl|/. The best fit values to 



Eq. (37) are given in the fourth column. The t-extensions of the lattices in 
the z~direction are also reported. 
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L, L, t Es Ef'^'^y 



4 


16 


20 


188.57(22) 




0.40 


5 


20 


30 


348.49(22) 




0.26 


6 


24 


30,40 


521.73(16) 




0.18 


7 


28 


40 


720.90(28) 




0.13 


8 


8 


30 


214.19(8) 




0.07 


8 


32 


30,40 


949.58(18) 




0.08 


9 


9 


30 


278.85(7) 




0.05 


9 


16 


20 


518.58(17) 




0.04 


9 


36 


30,40 


1207.47(17) 




0.08 


9 


64 


30 


2173.99(53) 




0.16 


10 


10 


30 


350.96(7) 


351.06 


0.04 


10 


40 


30,40 


1496.65(19) 


1496.47 


0.06 


11 


11 


30 


430.66(8) 


430.62 


0.03 


11 


44 


30,40 


1816.17(25) 


1815.98 


0.05 


12 


12 


30,40 


517.94(9) 


517.93 


0.03 


12 


48 


40 


2165.83(35) 


2166.18 


0.05 


13 


13 


30 


613.07(9) 


612.95 


0.02 


14 


14 


30,40 


715.73(17) 


715.66 


0.02 


15 


15 


30 


825.83(19) 


826.04 


0.02 


16 


16 


30,40 


944.16(18) 


944.08 


0.02 


18 


18 


30,40 


1203.11(15) 


1203.12 


0.01 


20 


20 


30,40 


1492.86(20) 


1492.73 


0.01 


22 


22 


30,40 


1812.86(21) 


1812.89 


0.009 


24 


24 


30,40 


2163.25(23) 


2163.58 


0.007 



Tab. 6. Data for the surface energy at (3 = 0.240 obtained by MC 
simulations with the local demon algorithm. The CWM best fit values of 



Eq. (31) to these data are given in the case Li > 10. In the last column the 



values of the two-loop parameter x defined in Eq. (|gqj are given. 



30 



/? 



approx. x^ /dof C.L. do/ 



a 



\/^ 



0.2240 


21 


0.63 


86% 


16 


0.004778(14) 


1.343(13) 




11 


2.02 


1 % 


16 


0.004839(14) 


1.561(15) 




cl 


66.0 


0% 


16 


0.004534(14) 


1.480(13) 


0.2246 


21 


0.54 


82% 


8 


0.006547(69) 


1.354(16) 




11 


1.89 


6% 


8 


0.006792(69) 


1.681(18) 




cl 


17.5 


0% 


8 


0.005987(73) 


1.539(16) 


0.2258 


21 


0.45 


89% 


8 


0.009418(61) 


1.271(14) 




11 


1.50 


15% 


8 


0.009587(61) 


1.511(15) 




cl 


20.2 


0% 


8 


0.008363(63) 


1.324(13) 


0.2275 


21 


1.05 


40% 


16 


0.014728(40) 


1.332(5) 




11 


2.36 


0% 


16 


0.015283(40) 


1.612(6) 




cl 


18.3 


0% 


16 


0.015114(39) 


1.594(6) 



Tab. 7. For each j3, we report the best fit results obtained using Eq. (37), 

Eq. f |7^ and Eq. f|7^, which correspond to 2-loop, 1-loop and classical 

approximation, respectively. 



scaling law approx. X /dof C.L. dof a^ 



\/V^ 



Eq. (H) 21 
11 

Eq. (^9|) 21 
11 


0.86 
0.94 
0.88 
0.94 


74% 44 1.47(1) 0.629(1) 1.333(8) 
56% 33 1.55(2) 0.633(1) 1.535(10) 
70% 44 1.32(1) 0.618(1) 1.331(8) 
56% 33 1.38(1) 0.622(1) 1.533(10) 


Tab. 8. Fits on 


the samples 


defined by Eq. 1^), with Xcut ~ 0.2 and 



Xcut ^ 0.15 for the two-loop Eq. ^68i) and the one-loop Eq. ^61[) 

respectively. 
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x'/dof 



scaling law 



Xn 



approx. 



C.L. dof Uc 



\/^ 



Eq. © 



Eq. (H) 



0.04 



0.07 



0.10 



0.12 



0.14 



0.16 



0.21 



0.04 



0.07 



0.10 



0.12 



0.14 



0.16 



0.21 



21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 
21 
11 



1.39 
1.50 
0.65 
0.72 
0.24 
0.35 
0.78 
0.82 
1.31 
1.96 
0.78 
0.74 
2.87 
6.16 
1.39 
1.50 
0.66 
0.72 
0.24 
0.34 
0.87 
0.88 
1.07 
1.70 
0.77 
0.73 
2.83 
6.19 



25% 
22% 



70 

61% 
95% 



70 

54% 
51% 
25% 
7% 



70 

62% 
1% 



70 

65% 
61% 



70 

48% 
47% 



70 

12% 



70 

62% 
1% 



2 
2 
5 
5 
5 
5 
4 
4 
6 
6 
6 
6 
7 
7 
2 
2 
5 
5 
5 
5 
4 
4 
6 
6 
6 
6 
7 
7 



1.53(6 
1.55(6 
1.49(6 
1.48(6 
1.48(4 
1.50(5 
1.44(6 
1.47(6 
1.49(4 
1.51(4 
1.50(7; 
1.50(7 
1.28(4 
1.58(5 
1.37(5 
1.39(5 
1.33(5 
1.33(5 
1.33(4 
1.34(5 
1.29(5 
1.32(5 
1.34(4 
1.35(4 
1.35(6 
1.35(6 
1.14(3 
1.40(4 



) 0.636(7; 


1.24(13) 


) 0.637(7; 


1.33(14) 


) 0.630(3; 


1.35(7) 


) 0.630(3; 


1.43(7) 


) 0.629(2; 


1.37(5) 


) 0.631(2; 


1.48(5) 


) 0.627(4; 


1.32(3) 


) 0.629(4; 


1.49(3) 


) 0.630(3; 


1.33(2) 


) 0.631(3; 


1.52(2) 


) 0.630(5; 


1.35(2) 


) 0.630(5; 


1.57(2) 


) 0.612(3; 


1.31(1) 


) 0.634(4; 


1.61(1) 


) 0.625(7; 


1.24(13) 


) 0.626(7; 


1.33(14) 


) 0.619(3; 


1.35(7) 


) 0.619(3; 


1.43(7) 


) 0.618(2; 


1.37(5) 


) 0.620(2; 


1.48(5) 


) 0.616(4; 


1.32(3) 


) 0.618(4; 


1.49(3) 


) 0.620(3; 


1.33(2) 


) 0.621(3; 


1.51(2) 


) 0.620(5; 


1.35(2) 


) 0.620(5; 


1.57(2) 


) 0.600(3; 


1.31(1) 


) 0.622(4; 


1.61(1) 



Tab. 9. Fits on the samples defined by Eq. (71). 
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L 



AEs 



AE^^^^ 



8 


-25.62(30) 


21.754(159) 


10 


-2.47(29) 


13.923(102) 


12 


3.79(25) 


9.669(71) 


14 


5.17(45) 


7.103(52) 


16 


5.42(36) 


5.439(40) 


18 


4.36(32) 


4.297(31) 


20 


3.79(39) 


3.481(25) 


22 


3.31(46) 


2.877(21) 


24 


2.58(58) 


2.417(18) 



Tab. 10. Differences of surface energies at (3 = 0.240, obtained from 

Eq. ^74 ) with the data of Tab. 6 and the notation L = 2Li = L2/2. The 

theoretical predictions A£^^ are given by Eqs. (T\,7i,'7^). 



Li> xii/dof xVdof C.L. dof 



a 



a 



\'/\ 



8 


107.2 


8.02 


0.00 


17 


0.0620(7) 


3.81383(24) 


33.10(8) 


9 


83.3 


2.25 


0.00 


15 


0.0593(8) 


3.81321(25) 


33.12(9) 


10 


48.5 


0.96 


0.48 


11 


0.0569(12) 


3.81292(29) 


33.17(13) 


11 


24.1 


0.66 


0.74 


9 


0.0565(18) 


3.81261(36) 


33.10(17) 


12 


8.4 


0.57 


0.78 


7 


0.0596(37) 


3.81213(50) 


32.83(27) 


13 


1.7 


0.72 


0.60 


5 


0.057(11) 


3.8123(14) 


32.95(87) 



Tab. 11. Data of Tab. 6 are fitted using Eq. (fTqj. For comparison, the 
reduced x^ 's are reported of the fits on the same sample of data using 



Eq. (6S). 



Li> 



a 



a 



a'/{2a) 



Eq. (|7|Q 
10 
11 



0.0590(2) 
0.0569(12) 
0.0565(18) 



3.813(2) 
3.81292(29) 
3.81261(36) 



32.31(13) 

33.5(7) 

33.7(1.1) 



0.593(2) 
0.615(13) 
0.619(20) 



Tab. 12. The scaling ratio a' /{2(t) is evaluated according to Eq. ^Tdj), with 

V ~ 0.63 for three estimates of cr and a' . In the last column the 

corresponding estimates of v, using Eq. ^^), are given. 



33 



Figure Captions 



Fig. 1. Histogram of the magnetization for a typical MC ensemble at 
(5 = 0.2275, with lattice sizes Li = 20, L2 = 23 and t = 120. 

Fig. 2. MC data and best fit curves obtained with the two-loop formula 
( p^ j plotted for different values of /3 versus the classical interface area A = 
L1L2. From top to bottom, the first three lines correspond to jS = 0.2246, 
the next three to (3 = 0.2258 and the last two to (3 = 0.2275. The values are 
also given in Tabs. 2, 3 and 4, respectively. 

Fig. 3. The same as Fig. 2 but for (3 = 0.2240. 

Fig. 4. x^ /dof of the fits made with the one-loop Eq. ( |g^ j (empty cir- 
cles) and two-loop Eq. (J6^ (full circles) on the samples defined by Eq. ([7^. 



Fig. 5. Values of the constant \l \fo obtained fitting on the samples 
defined by Eq. ^7l\) with the one-loop Eq. ( fg^j (empty circles) and two-loop 
Eq. M^ (full circles). 



Fig. 6. Comparison between Monte Carlo data and theoretical predic- 
tions (see Eq. ( f7^ in the text) for the surface energy differences defined in 
Eq. (f^[j. The corresponding values are also given in Tab. 10. The dotted 
line is the theoretical expectation if two-loop corrections are neglected. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407002vl 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407002vl 
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This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407002vl 



fig. 3 
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This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407002vl 
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This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407002vl 
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This figure "figl-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407002vl 
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